Surface-edge state and half quantized Hall conductance in topological insulators 
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We investigate the transport properties of the surface states of a three-dimensional topological 
insulator in the presence of a spin-splitting Zeeman field. We propose a picture that the chiral edge 
state forms on the surface, and is split into two halves that are spatially separated, each carrying one 
half of the conductance quantum (e 2 /h). This picture is confirmed by numerical simulation in a four- 
terminal setup. It is demonstrated that the difference between the clockwise and counterclockwise 
transmission coefficients of the two neighboring terminals is approximately one half, which suggests 
that the half quantized Hall conductance can be manifested in an appropriate experimental setup. 
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Topological insulators are insulating in the bulk, but 
have metallic surface states possessing an odd number of 
Dirac cones of massless fermions pj[3] . The band struc- 
ture and the quantum spin texture of these surfaces states 
have been well established theoretically and experimen- 
tally [IHS]. In the presence of a spin-splitting Zeeman 
field, which could be induced by magnetically doping the 
samples or putting the samples in the proximity of ferro- 
magnetic materials, the Dirac fermions will gain a mass 
and the spectrum opens a gap [101 E] • When the Fermi 
level is located within the energy gap, it was proposed 
that the Hall conductance of the surface states will be 
one half of the conductance quantum e 2 /h [T2HI4] . Based 
on this, Qi et al. proposed the unconventional magneto- 
electric effect, which is regarded as one of the character- 
istic features of the topological insulators [T4HT6] . 

On the other hand, it is not clear whether or not the 
half quantization of the Hall conductance can be directly 
observed in the transport measurement. In the usual 
quantum Hall system, the current-carrying chiral edge 
states are responsible for the integer quantized Hall con- 
ductance measured in the transport experiment. [T71 [T8] 
It is not immediately clear whether or not the similar chi- 
ral edge state will form on the closed surface of a topolog- 
ical insulator, and how the quantized nature of the edge 
states can be reconciled with the prediction of the half 
quantization of the Hall conductance. [T2HT4] 

To get a definite answer to these questions, we inves- 
tigate the multi-terminal transport properties of a 3D 
topological insulator in the presence of a uniform spin- 
splitting Zeeman field. We propose that the the closed 
surface of the topological insulator will form two insu- 
lating domains of the different topologies (i.e., positive 
vs. negative gaps), separated by a gapless metallic belt. 
A chiral edge state will form and is concentrated around 
the boundaries between the insulating domains and the 
metallic belt. Effectively, the chiral edge state is split 
into two halves, each of which is circulating around the 
boundary of one of the domains and carrying one half of 



the conductance quantum (e 2 /h). Such a picture recon- 
ciles the apparent conflict between the half quantization 
and the index theorem. It also suggests that by attaching 
the electrodes to the boundary of one of these domains, 
it is possible to directly measure the half quantized Hall 
conductance. The picture as well as its prediction are ver- 
ified by the numerical simulation based non-equilibrium 
Green's function method. 
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FIG. 1: (Color online) (a) Schematic of a 3D TI in a weak 
Zeeman field, and the formation of chiral current on the top 
and bottom surface boundaries, (b) A chiral edge state will 
form around the domain wall between the 2D Dirac fermions 
with positive and negative masses, and the wave function is 
illustrated. The arrow indicates the flow of edge current, (c) 
When the sharp domain wall evolves to finite-width metallic 
band, the edge mode is effectively split to two halves concen- 
trated around the two boundaries. 

To illustrate the basic physics, we consider a 3D topo- 
logical insulator of the cubic shape. A Zeeman field is 
applied along the z-direction, as shown in Fig. 1(a). Be- 
cause the bulk of the system is insulating, it is effectively 
a closed 2D surface with six facets. The effective Hamil- 
tonian of the Dirac fermions for the surface state can be 
written as [T91I20] : 

H e fi(k) =v(kx a)-n-g\\fi B h\\a\\ -g±fiBh±-a±, (1) 



where n denotes the normal vector of the surface, a = 
{cr X: (7y : a z } are the Pauli matrices, h\\ (cry) and hj_ (crj_) 
are the Zeeman field (Pauli matrix) components paral- 
lel and perpendicular to the normal vector, respectively, 
and #|| and g± are the corresponding spin ^-factor. Note 
that the surface state has the anisotropic spin ^-factor 
due to the strong spin-orbit coupling of the bulk band: 
g\\ is the same as the ^-factor of the bulk material, and 
g± is renormalized by bulk band parameters and is usu- 
ally strongly suppressed [19] [20]. The different facets 
of the surface have the different effective Hamiltonians 
respective to the different normal vectors n. For the 
top and the bottom facets, the effective Hamiltonian can 
be written as H e & = ±v(k x a y — k y a x ) + A z a z , where 
+ and — are for the top and bottom surfaces, respec- 
tively, and A z = —gnfish. The spectrum will open 
a gap on these facets, and the Dirac fermions gain a 
mass ±A Z . On the other hand, the effective Hamilto- 
nians for the side facets can be written as the forms like 
#eff = v[(k x + Ak)a z - k x a x ], where Ak = g^fish. In 
this case, the Zeeman field has no effect except shifting 
the Dirac point. When the fermi level is located in the 
gap of the top and bottom surface, the system becomes 
effectively two insulating domains separated by a con- 
ducting belt with massless Dirac fermions. 

The two insulating domains belong to the different 
topological classes due to the fact that the normal vectors 
for the top and the bottom surfaces are in the opposite 
directions. We consider the limiting case that the sample 
has the vanishing thickness along the z-direction (but ig- 
nore the coupling between the top and bottom surfaces). 
The system is locally equivalent to a domain-wall struc- 
ture across which the Dirac fermion masses changes from 
positive to negative, as shown in Fig. 1(b). The change 
of the band topology necessitates the close of the energy 
gap around the domain wall, giving rise to the gapless 
edge mode. By directly solving the wave equation and 
matching the wave functions at the boundary, it is easy 
to obtain that the wave function of the edge mode has the 
form ^(a;,*/) = V / A^/4^(l, 1) T exp[-(A z /v)\x\ + ik y y], 
and the dispersion is linear in the momentum: E = 
—vk y . The edge mode will form an ideal one-dimensional 
chiral edge channel carrying a quantized conductance 
e 2 /h. The channel cannot be backscattered and is robust 
against the disorder. Transport measurement will yield 
a quantized Hall conductance if attaching electrodes to 
the edge of the sample. 

When the thickness of the sample along the ^-direction 
is finite, the two insulating domains are separated by a 
finite metallic region of massless Dirac fermion, equiv- 
alent to the configuration shown in Fig. 1(c). In this 
case, the electrons in the metallic region will in gen- 
eral form sub-bands due to the confinement by the in- 
sulating domains. The normal sub-band has dispersion 



momentums of the eigen- modes along the confinement di- 
rection. Beside these, there is a chiral solution with the 
dispersion E = —vk yi corresponding to k n = but with 
only the left-going mode. The chirality is also shown in 
the penetration amplitude of the normal sub-band wave 
function into the insulating domains: the penetration 
amplitude maximizes when E ~ —vk y: and decreases 
when moving away from it. All sub-bands crossing the 
Fermi energy will contribute to the transport properties, 
and form a novel chiral edge state as a result of collective 
combination (not a single mode). The chirality will give 
rise to a finite Hall conductance. 

To see what the transport measurement will yield for 
such a system, we employ the direct numerical simula- 
tion. We investigate an isotropic 3D topological insulator 
with the tight-binding Hamiltonian on a cubic lattice: 

H = ^c\MoCi+ Yl (4TaCi + a + 4 +a T a C^ (2) 



z,a=x,y,x 



where T a = Ba z ® <Jo — i^cr x <g> a a and Mq = [M — 
6B)<j z <S) (Jo + A 2 <jo (8) cr z . The lattice space is taken to 
be unity. Near the Y point in the k space, it is reduced 
to a Dirac-like model, which can be derived from either 
the theory of invariants [7 , or the 8-band extended Kane 
model near the inversion of the Tq and Tg" bands. The 
Dirac model for a topological insulator has been discussed 
by several authors [HI [21] . It was shown that the model 
yield the effective surface Hamiltonian Eq. [I] with v = 
A p~9j [20] . We first employ the Green function technique 
to calculate the local density of states (LDOS) for the 
surface states, in particular the spatial distribution near 
the boundary, to demonstrate the formation of the edge 
state. Then we calculate the Hall conductance for a four- 
terminal setup, which directly reflects the output of the 
transport measurement. 
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E n = ±vJk 2 + &2 , where k n ^ is one of the discrete 
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FIG. 2: (Color online) Local density of state on an infinite xy 
surface of a semi-infinite 3D system, (left) gapless single Dirac 
cone of the surface state; (right) gap opening by application of 
a Zeeman splitting term. The model parameters are A — 0.5, 
B = 0.25, M = 0.3, and A z = 0.07. 

In Fig. 2 we present LDOS p(k x , k y ) of the top surface 
for a semi-infinite sample with z G [0, — oo). In this case, 



k x and k y are good quantum numbers. A recursive ap- 
proach are employed to calculate the Green's function for 
a few top layers [22 . LDOS at the top surface can be 
calculated by pi(k x ,k y ) = — -Ti[ImGu(k x ,k y )], where 
Gu is the retarded green function for the ith. layer from 
the top. In the absence of the Zeeman field, the gapless 
linear Dirac cone can be observed. Position dependence 
of LDOS along the z-axis shows that the states in the 
Dirac cone reside dominant ly near the top surface, which 
is characteristic of the surface states for a topological in- 
sulator, and is consistent with the analytical solution. [19] 
In the presence of the Zeeman field, the LDOS shows that 
a gap of magnitude 2A 2 is opened, as expected. 
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FIG. 3: (Color online) LDOS on the top surface of a structure 
that is infinite in X, finite in Y and semi-infinite in Z direction. 
Sampling is taken correspondingly in a, b and c regions as 
illustrated in the upper panel. A z = 0.15, M = 0.4, and L y = 
30a (a is the lattice space). 

To see the distribution of the surface states near the 
edge of the top surface, we study a geometry that is fi- 
nite in the y axis, infinite in the x axis and semi-infinite 
in the z axis, as illustrated in the upper panel of Fig. 
3. In this case only k x is a good quantum number. The 
LDOS is p(y,k x ) = -^Tr[ImG 00 (y,k x )}. The average 
LDOS at three small different regions (a)-(c) along the 
^-direction are shown in the lower panel of Fig. 3. It can 
be seen that the finite density of states emerges in the 
gap at the locations (a) and (c) near the boundaries of 
the top surface. The finite density of the in-gap states 
originates from the penetrated wave functions from the 
side facets, which will decay exponentially in space when 
the energy of the state is located in the gap of the top 
surface states. The chiral nature of these in-gap states 
can be directly observed: the LDOS at the edge (a) max- 
imizes at E = — vk x (left going), while the LDOS at the 
edge (c) maximizes at E = vk x (right going). As a result, 



these states will carry a net chiral edge current. 

We calculate the current density profile of these in-gap 
states at the fermi surface Ef along the ^/-direction us- 
ing (J x (y)} = ie f kx Tr[v x (r,k x )G < (r,k x )} where v x = 
dH/dk x is the velocity operator and Ef is in the Zee- 
man splitting gap. In an equilibrium condition we have 
G< = f(E f )(G a - G r ), where f(E f ) is the Fermi-Dirac 
distribution function and we assume zero temperature. 
In Fig. 4a the current density distribution along the y- 
direction for the top five layers is shown. It can be seen 
that the current is localized dominantly near the edges 
in the insulating top layers. To see the distribution of 
the chiral edge current in the side facets, we sum up the 
current density for each layer over the half the width of y: 
J t = J^q y J x (y) and observe its layer dependence. As is 
shown in Fig. 4b, J t shows the typical behavior of Friedel 
oscillation [23] , which maximize at the boundary and de- 
cays to zero when moving away from the boundary. We 
see that the distribution of the chiral edge current is con- 
centrated around the boundary, even the wave- functions 
in the metallic region are all extended. 



0.6 
0.3 
_* 
-0.3 
-0.6 
0.3 
0.0 
-t; -0.3 
-0.6 
-0.9 



r-J 


1 1 ■ 1 ■ 

: 1 


i 1 i 1 i 


i i i i i 


1 1 1 1 1 


: 


:"" ■_ 


r j 

i_ _ 


- 0) - 


",1,1," 


' " 


" 


" 





10 20 





10 20 


10 20 





10 20 





10 20 


' i 

i 




i 


1 i ' 


1 






1 i 

(b) - 




m mi g rg i n m » 



0.2 0.25 0.3 0.35 
Ef2 



FIG. 4: (Color online) (a) Local current density J x distribu- 
tion along the width of Y (L y = 30a) for the top 5 layers 
(the leftest is the 1st); (b) Total chiral current density J t as 
a function of layer depth Z; (c) Transmission coefficients of 
the 4 terminal device, Ef 2 is fixed, the dashed line indicates 
the gap position; (d) same with (c), Ef\ is fixed, the insert 
schematically illustrated the device set up. Other parameters 
are A z = 0.15, M = 0.4, E f = 0.075. 



After establishing the existence of the chiral surface- 
edge states, we calculate the Hall conductance numer- 
ically using the Landauer-Biittiker formalism [241126] . 
The set up of the device is illustrated in the insert of 
the Fig. 4d: four 2D metallic leads (/i = 1,2,3,4) are 
attached to the top surface of a semi-infinite 3D topo- 
logical insulator, acting as the measurement electrodes. 
The Zeeman field is normal to the top surface. The multi- 
terminal conductance can be deduced from the transmis- 
sion coefficient T pq from the terminal p to terminal g, 
T pq = Tr[T p G R r q r q ] where T p is determined by the self 
energy at the terminal p [25 . The calculated transmis- 
sion coefficients as a function of chemical potentials E\f 
in the 3D topological insulator and E 2 f in the terminals 
are plotted in Figs. 4c and 4d. When the chemical po- 
tential Eif is located in the gap of the surface states, the 
transmission coefficients show the "half chirality", i.e., 
T pq — T qp ^ ±1/2 between the two neighboring termi- 
nals p and g, and T pq = between the non-neighboring 
terminals. This is different from the chirality shown in 
the usual quantum Hall system, where T pq — T qp = ±1 
between the two neighboring terminals. It is also impor- 
tant to note that the transmission coefficient between the 
two neighboring terminals could be very large due to the 
presence of the metallic side facets, and it is the differ- 
ence between the clockwise and anti-clockwise transmis- 
sion coefficients that shows the " half quantization" . 

The predicted half quantization can be directly mea- 
sured. The straightforward way to observe it in the four- 
terminal setup is to apply a voltage between the termi- 
nals 1 and 3 (V13), and measure the current between the 
terminal 2 and 4 (h^)- It is easy to show that the cross- 
conductance cr 2 4,i3 = WVia = (e 2 /2h)(T 12 - T21), 
yielding e 2 /Ah for the half quantization. The measure- 
ment using the usual six-terminal Hall bar configuration 
could be more tricky due to the presence of the metal- 
lic side facets, which give rise to the finite longitudi- 
nal conductance cjl. In the limit of thick sample with 
o~l ^> e 2 //i, the Hall conductance an should approach to 
(4e 2 /h)(T 12 - T 2 i). It yields 2e 2 /h for the half quanti- 
zation. It can be distinguished from the usual quantized 
Hall effect by the finite longitudinal conductance. 

To summarize, we find a chiral surface-edge state for 
a 3D topological insulator in the presence of Zeeman 
splitting. Effectively, one can consider that the origi- 
nal chiral edge mode carrying one conductance quantum 
e 2 /his split into two halves concentrated around the two 
boundaries between the insulating domains and metallic 
region. Its origin is quite different from those in inte- 
ger quantum Hall and quantum spin Hall effect, in which 
there exist well defined chiral edge states. Here, many 
conducting channels, each contributing a fraction of the 
chirality, may give rise to an effective edge state concen- 
trated around the domain boundaries, as indicated by 
the spatial distribution of the chiral edge current. The 
half quantization can be measured in the multi-terminal 



setup. Moreover, the existence of the effective edge state 
may be directly observable. For example, because the 
spin polarization of Dirac fermions in the surface states 
is proportional to the charge current, the spatial oscilla- 
tory pattern of the chiral edge current associating with 
the effective edge state could be directly probed by the 
spatially resolved Kerr rotation technique [27] . 
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